Accuracy of the pion elastic form factor extracted from a local-duality sum rule 
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We analyze the accuracy of the pion elastic form factor predicted by a local-duality (LD) version of 
dispersive sum rules. To probe the precision of this theoretical approach, we adopt potential models 
with interactions that involve both Coulomb and confining terms. In this case, the exact form factor 
may be obtained from the solution of the Schrodinger equation and confronted with the LD sum-rule 
results. We use parameter values appropriate for hadron physics and observe that, independently of 
the details of the confining interaction, the deviation of the LD form factor from the exact form factor 
culminates in the region Q 2 « 4-6 GeV 2 . For larger Q 2 , the accuracy of the LD description increases 
rather fast with Q 2 . A similar picture is expected for QCD. For the pion form factor, existing data 
suggest that the LD limit may be reached already at the relatively low values Q 2 — 4-10 GeV 2 . Thus, 
large deviations of the pion form factor from the behaviour predicted by LD QCD sum rules for higher 
values of Q 2 , as found by some recent analyses, appear to us quite improbable. New accurate data on 
the pion form factor at Q 2 = 4-10 GeV 2 expected soon from JLab will have important implications 
for the behaviour of the pion form factor in a broad Q 2 range up to asymptotically large values of Q 2 . 

PACS numbers: 11.55.Hx, 12.38.Lg, 03.65. Ge, 14.40.Be 
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1. INTRODUCTION 

In spite of a rather long history of theoretical studies of the pion elastic form factor, no consensus on its behaviour in 
the spacelike region Q 2 > 2-4 GeV 2 has been reached so far. For instance, recent theoretical investigations [IHH report 
results for the pion form factor much larger than earlier ones (cf. Fig. [J): According to [H-i], even at Q 2 « 50-100 GeV 2 
the pion form factor F n (Q 2 ) remains much larger than the asymptotic behaviour expected from perturbative QCD Q 



Q 2 F 7r (Q 2 ) = SirasiQ 2 )/, 



(1.1) 



with the decay constant of the pion and a s (Q 2 ) the strong coupling. Subleading logarithmic and power corrections 
modify the behaviour (jl.ip at large but finite Q 2 . In early applications of QCD one hoped that power corrections would 
vanish fast enough with Q 2 ; however, later investigations revealed that nonperturbative power corrections dominate the 
form factor F 7T (Q 2 ) up to relatively high Q 2 w 10-20 GeV 2 . This picture has arisen from different approaches 
At this stage, the conclusion was that even for Q 2 as large as Q 2 — 20 GeV 2 the O(l) term provides about half of the 
form factor and the perturbative-QCD formula based on factorization starts to work well only at Q 2 > 50-100 GeV 2 . 
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Fig. 1: Some predictions for the pion elastic form factor F n {Q 2 ) — lower solid (black) line: BLM'2008 0], upper solid (green) line: 
BT'2008 [J], dashed (magenta) line: GR'2008 ;2], and dash-dotted (blue) line: BPS'2009 5, Eq. (4.11b)] — vs. experiment 
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A convincing and largely model-independent argument comes from QCD sum rules in their local-duality (LD) version 
(l5j . which enables one to take into account, on an equal footing, the 0(1) and the 0(a s ) contributions to F^(Q 2 ): the 
LD form of QCD sum rules predicts in an essentially model-independent way the relative weights of these contributions 
|H, but needs an additional input — the <5 2 -dependent effective threshold s e g(Q 2 ) — in order to calculate F n (Q 2 ). All 
previous applications of LD sum rules had to rely on assumptions about the Q 2 behaviour of the effective threshold. 
The LD model [l5|,[l6| assumes that reasonable estimates for the form factor already starting from Q 2 > 2-4 GeV 2 may 
be obtained by setting s c s (Q 2 ) equal to its LD value sld = 47r 2 / 2 . Under this assumption, Ft T (Q 2 ) has been calculated 
in i. 1 

Surprisingly, recent studies [11-0, H| find large deviations from the LD results Q. To quantify these deviations, Fig. [2] 
compares the equivalent effective thresholds 2 computed from the results of these analyses. For some of these studies, 
the deviation of the equivalent threshold from sld rises with Q 2 , although the accuracy of the LD model is expected 
to increase with Q 2 . However, all of these studies involve explicit or implicit assumptions. 

s ef f(Q 2 ) [GeV 2 ] 




Fig. 2: Equivalent effective thresholds s e g(Q 2 ) resulting from the findings of P,0,I3|, identified by the same line codes as in Fig.[T] 

The goal of this paper is to test the accuracy of the LD model by exploiting the fact that in quantum mechanics the 
form factor may be found in two ways, viz., by the LD sum rule and by solving the Schrodinger equation. For the first 
time we calculate the exact <5 2 -dependent threshold in a quantum-mechanical potential model for different confining 
interactions. We find that the accuracy of the LD approximation at a given value of Q 2 depends on the details of the 
confining interaction. More importantly, we observe a universal feature that does not depend on these details: namely, 
for realistic values of the parameters of the potential model relevant for hadron physics, the accuracy of the LD model 
increases with Q 2 already starting at relatively low values Q 2 « 4-6 GeV 2 . 

In QCD, experimental data on the pion form factor at small Q 2 indicate that the LD limit may be reached already at 
relatively low values of Q 2 ~ 5-6 GeV 2 . This, of course, does not mean that the asymptotic pQCD formula describes 
well the form factor in this region: the O(l) term still dominates the form factor up to Q 2 « 20 GeV 2 . Therefore, we 
conclude that large and growing-with-Q 2 deviations from the LD limit in the region Q 2 m 20-50 GeV 2 implied by the 
recent analyses [l|-|3| seem to us very unlikely. 

This paper is organized as follows: In the next section, after recalling rather well-known basics of the LD limit of the 
QCD sum-rule approach we formulate, for use in these LD sum rules, a model for s e ff (Q 2 ) slightly more sophisticated 
than the one used in Q. In Sec. [31 we discuss in detail rigorous features and assumptions employed in the LD approach. 
Section [4] presents the quantum-mechanical counterpart of the LD model for F 7r (Q 2 ) using potentials that involve both 
Coulombic and confining interactions: we examine the impact of the specific shape of the confining interactions on the 
accuracy of the LD model for the elastic form factor. Section [3] summarizes our conclusions and outlook. Appendix |A1 
presents several technical details of perturbative two-loop calculations in quantum mechanics. 



1 Notice that the pioneering work on the LD model [lf| used a simple interpolation formula for the unknown 0(a s ) contribution to -FV(Q 2 ). 
The first analysis of the pion form factor taking into account the 0(a B ) contribution obtained in [l7ll was performed in [5|. 

2 The equivalent effective threshold is defined be requiring that a given theoretical prediction for the pion form factor is reproduced by the 
LD sum rule 112.91 1 if the corresponding equivalent threshold is used (cf. Sect. 2). 
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2. SUM RULE 

The fundamental objects for a sum-rule extraction of pion features are the two- and three-point correlation functions 

n( P 2 ) = J d 4 xe^ (n \T[j(x) jt( )]| n> , 

T(Pi,pW) = [ d 4 x 1 d 4 x 2 e i ^'^{n\T[j(x 1 )J(0)jHx 2 )]\n), q = Pi-p 2 , Q 2 = -<? 2 ; (2.1) 



where f2 labels the physical vacuum, j is shorthand for the interpolating axial current j§ a of the positively charged pion, 
|j5a(0)| 7r(p)) = ipafn, J labels the electromagnetic current J^, and for brevity we omit all Lorentz indices. In QCD, 
the correlators (I2.1j) can be found by applying their OPEs. Instead of discussing the Green functions (|2.ip in Minkowski 
space, it is convenient to study the time-evolution operators in Euclidean space, which arise upon performing the Borcl 
transformation p 2 r to a parameter r related to Euclidean time. The Borel image of the two-point correlator II (p 2 ) is 

oo 

IIope(t) = J dse _ST p pee t(s) +n cond (r), p pcr t(s) = p (s) + a s pi(s) + 0{a 2 s ), (2.2) 
o 

with spectral densities Pi(s) related to perturbative two-point graphs, and nonperturbative power corrections n con d(r). 
At hadron level, insertion of intermediate hadron states casts the Borel-transformed two-point correlator into the form 

n(r) = fl e-' n * T + excited states. (2.3) 

In this expression for II(t), the first term on the right-hand side constitutes the pion contribution. Applying the double 
Borel transform p\ 2 — > r/2 to the three-point correlator T{p\ 1 p\ 1 q 2 ) results, at QCD level, in 

oo oo 

Tope(t, Q 2 ) = J J dsi ds 2 exp (~ Sl ^ 52 r^j A pclt (s 1 ,s 2 ,Q 2 ) + r cond (r, Q 2 ), 
o o 

A P crt(si, 32, Q 2 ) = A (si, s 2 , Q 2 ) + a s Ai(si, as, Q 2 ) + 0(q 2 ), (2.4) 

where A pert (si, s 2 , Q 2 ) is the double spectral density of the three-point graphs of perturbation theory and r conc i(T, Q 2 ) 
labels the power corrections. Inserting intermediate hadron states yields, for the hadron-level expression for T(t, Q 2 ), 

T(t, Q 2 ) = F 7I (Q 2 ) f 2 e" m ' T + excited states. (2.5) 

Quark-hadron duality assumes that above effective continuum thresholds s c g the excited-state contributions are dual 
to the high-energy regions of the perturbative graphs. In this case, the relevant sum rules read in the chiral limit [l5l[20j 



Seff(T) 



(a s G 2 ) 176tt a s (qq} 2 
12tt T+ 81 



f 2 = ds e— p pc As) + r + — ^ r 2 + ■ ■ • , (2.6) 



o 



F n (Q 2 )f%= j j d Sl d S2 Ap Crt ( S i, S2 ,Q 2 )exp^-^ 



+ <^ r+ 4 ! r^ r2(i3 + g2r)+ __ (2?) 

As a consequence of the use of local condensates, the right-hand side of (|2.7p involves polynomials in Q 2 and therefore 
increases with Q 2 , whereas the form factor F 7T (Q 2 ) on the left-hand side should decrease with Q 2 . Therefore, at large Q 2 
the sum rule (|2.7I) . with its truncated series of power corrections, cannot be directly used. There are essentially only two 
ways for considering the region of large Q 2 . 

One remedy is the resummation of all power corrections: the resummed power corrections decrease with increasing 
Q 2 . This may be achieved by the introduction of nonlocal condensates [l8[ in a, however, model-dependent manner [lj| . 

Another — rather simple — option is to fix the Borel parameter r to the value r = 0, thus arriving at a local-duality 
sum rule [l5j . Therein all power corrections vanish and the remaining perturbative term decreases with Q 2 . In the LD 
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limit, one finds 



Seff 



fl = J ds Ppcit ( S ) = |2| (l + ^) + 0(al) , (2.8) 
o 

s ff(Q 2 ) s cff (Q 2 ) 

F^{Q 2 )fl= J J d Sl d S2 A port ( Sl , S2 ,Q 2 ). (2.9) 



The spectral densities p per t(s) and A pcrt (s 1 , s 2 , Q 2 ) are calculable by perturbation theory. Hence, by fixing s e s and 
s c s(Q 2 ), it is straightforward to extract the pion's decay constant f„ and form factor F 1T {Q 2 ). 

Noteworthy, the effective and the physical thresholds are different quantities: The latter is a constant determined by 
the masses of the hadron states. The effective thresholds s c g and s e e are parameters of the sum-rule method related to 
the specific realization of quark-hadron duality; in general, they are not constant but depend on external kinematical 
variables pi |22T| . 

Let us recall the important properties of the spectral densities on the right-hand sides of (|2.2p and (|2.4p : For Q 2 — » 0, 
the Ward identity relates the spectral densities Pi(s) and Aj(si, s 2 , Q 2 ) of two- and three-point functions to each other: 

lim A i (s 1) s 2 ,Q 2 ) =p i (s 1 )S(s 1 -s 2 ), i = 0,l,.... (2.10) 

Q 2 ->o 

For Q 2 —> oo and Sx,2 kept fixed, explicit calculations [l?} yield 

1 &7T 

lim A (s 1 ,s 2 ,Q 2 ) oc — j, lim A 1 (s 1 , s 2 , Q 2 ) = — -s Po(si) PoM- (2-11) 

Q 2 —¥oo Q 2 — >oo 

For the pion form factor Ft T (Q 2 ) on the left-hand side of (|2.7[) two exact features are known, namely, its normalization 
condition related to current conservation, requiring F^(0) = 1, and the factorization theorem (jl.ip . Obviously, if we set 

47T 2 f2 

Soft (Q 2 -> 0) = — , Seff (Q 2 ^ Co) = SLD = 47T 2 / 2 , (2.12) 

1 + a s (0)/7r 

the form factor F 7T (Q 2 ) extracted from the LD sum rule (|2.9[) satisfies both of these rigorous constraints. At small Q 2 , 
we assume a freezing of a s (Q 2 ) at the level 0.3, as is frequently done. 

With these preliminaries at our disposal, we are now in the position to formulate our approach to the pion elastic form 
factor Fir(Q 2 ) within the framework of QCD sum rules in the LD limit: 

1. We extract F 7r (Q 2 ) from the dispersive three-point QCD sum rule in LD limit (12.71) : the latter has the favourable 
feature that all power corrections vanish and all details of nonperturbative dynamics are encoded in one quantity, 
the effective threshold s e s (Q 2 )- We take into account the perturbative spectral densities up to 0(a s ) accuracy. 

2. It is easy to construct a model for s c g(Q 2 ) by some smooth interpolation between its values at Q 2 = 0, defined by 
the Ward identity, and Q 2 — > oo, determined by factorization: a simple parametrization with a single constant Qq 
fixed by fitting the data at Q 2 = 1 GeV 2 might read 



Soff(<9 2 



47T 2 / 2 



1 + tanh 



(0) 



Q 2 = 2.02 GeV^. (2.13) 



l + a s (0)/7r 

According to Fig.[3j our interpolation perfectly describes the well-measured data in the range Q 2 w 0.5-2.5 GeV 2 . 

Note that the effective continuum threshold s e g (Q 2 ) in Eq. (|2 . 1 3[) approaches its limit sld already at Q 2 sa 4-5 GeV 2 . 
For Q 2 > 4-5 GeV 2 , it practically coincides with the LD effective threshold of [l5[ . Moreover, for Q 2 > 5-6 GeV 2 the 
formula (|2.13|) is pretty close to the model of [4(. 3 Hence, also the resulting prediction for F V (Q 2 ) is rather close to the 
one we found earlier [4] (cf. Fig. [TJ. Obviously, the model labeled BLM in Fig. [1] provides a perfect description of the 
available F 7T (Q 2 ) data in the region Q 2 = 1-2.5 GeV 2 . For Q 2 > 3-4 GeV 2 , it reproduces well all the data, except for a 
point at Q 2 — 10 GeV 2 , where it is off the present experimental value, which anyhow has a rather large error, by some 
two standard deviations. 4 



3 It goes without saying that our present goal is not to improve the model of Q but to probe the accuracy of this model as a function of Q 2 . 

4 It is virtually inconceivable to construct models compatible with all experimental results within Q 2 = 2.5-10 GeV 2 , as revealed by closer 
inspection of Fig.[l] those approaches which hit the data at Q 2 = 10 GeV 2 overestimate the better-quality data points at Q 2 fa 2-4 GeV 2 . 
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Fig. 3: Effective continuum threshold s E s(Q 2 ) as function of Q 2 . Bold (red) line: exact behaviour s e a{Q 2 ) as reconstructed from 



the experimental data for (Q 



for Q 2 < 2.5 GeV . Solid (black) line: our simple interpolating model (|2. 13f) for s e ff (<2 



3. WHAT IS PREDICTED BY AND WHAT IS CONJECTURED IN THE LD MODEL? 



Interestingly, in the region Q 2 > 3-4 GeV 2 the BLM model yields considerably lower predictions than the results of 
the different theoretical approaches presented in Refs. [H-0,Q. For the F 7r (Q 2 ) predictions of 1, 2, 5], the corresponding 
equivalent effective thresholds s e s(Q 2 ) recalculated from (|2.9p are depicted in Fig. [2] In all instances, they considerably 
exceed, for larger Q 2 , the LD limit sld dictated by factorization. Moreover, their deviation from sld increases with Q 2 . 

Let us inspect more carefully what is, in fact, predicted by the LD sum rule and what is conjectured in this approach. 

The sum rule (|2.9[) for the pion form factor relies on two ingredients: first, on the rigorous calculation of the spectral 
densities of the perturbative-QCD diagrams (recall that power corrections vanish in the LD limit r = 0); second, on the 
assumption of quark-hadron duality, which claims that the contributions of the hadronic continuum states may be well 
described by the diagrams of perturbation theory above some effective threshold s c g . Thus, the only — although really 
essential — ingredient of the LD sum rule for the pion elastic from factor is this effective continuum threshold s c ff(Q 2 )- 
Let us emphasize that, since the 0(1) and 0(a s ) contributions to the pion form factor are governed by one and the same 
effective threshold s c s (Q 2 ), the relative weights of these contributions may be predicted. Their ratio Fi 0) (Q 2 )/^ (Q 2 ) 
turns out to be relatively stable with respect to s e ff and may therefore be calculated relatively accurately (see Fig. [U ■ 



F^ (0) (Q 2 , Seff) 
F„ (1) (Q 2 , Seff) 

3 r- 




Fig. 4: Ratio F^ 0) (Q 2 ) / F^ ] (Q 2 ) of O(l) and 0{a s ) contributions to the pion elastic form factor vs. Q 2 for two reasonable values 
of the — here by assumption constant — effective threshold s e ff. 
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Quark-hadron duality implies that the effective threshold (although being a function of Q 2 ) always — that is, also 
for large Q 2 — stays in a region close to 1 GeV 2 . Moreover, in order to satisfy the QCD factorization theorem for the 
0(a s ) contribution to the form factor, the effective threshold should behave like s g(Q 2 ) — > 47r 2 / 2 for Q 2 — > oo. This 
requirement has immediate consequences for the large-Q 2 behaviour of different contributions to the pion form factor: 

1. Since s c g(Q 2 ) is bounded from above, theO(l) contribution F^ \Q 2 ) to the elastic form factor F^(Q 2 ) of the pion 

behaves like (Q 2 ) oc 1/Q 4 for Q 2 — > oo. We would like to emphasize that the decrease of the soft contribution 
to the pion elastic form factor like l/Q 4 is a direct consequence of perturbation theory and quark-hadron duality. 

2. Consequently, for large Q 2 the pion elastic form factor F 7r (Q 2 ) is dominated by the 0(a s ) contribution F^ (Q 2 ). 

Now, in the holographic models of pl@] merely the soft contribution is considered: it behaves like l/<5 2 for large Q 2 . 
This is only possible if the effective threshold s c h(Q 2 ) rises with Q 2 . However, this immediately leads to the violation of 
the factorization theorem for the 0(a s ) contribution, which is governed by the same effective threshold. Consequently, 
we would like to emphasize that the findings of [l|, [|[ would imply that the QCD factorization theorem is violated. We 
thus conclude that the predictions of [3, HJ for the pion form factor seem to us improbable. 5 

Needless to say, the conventional LD model [l6| for the effective continuum threshold s e g (Q 2 ), defined by the choice 
Scff (Q 2 ) = 47r 2 / 2 for all Q 2 , or the slightly more sophisticated approach of [4| are approximations which do not account 
for all the subtle details of the confinement dynamics. Consequently, it is of utmost importance to acquire a satisfactory 
understanding of the accuracy to be expected within this approach, in other words, to obtain some reliable estimate of 
the expected deviations of the exact s e g(Q 2 ) from its LD limit sld in the momentum region Q 2 > 4-6 GeV 2 . To this 
end, let us take advantage of the fact that in quantum mechanics all the bound-state features may be found exactly by 
solving the Schrodinger equation. On the other hand, also within quantum mechanics we may construct LD sum rules. 



4. LOCAL-DUALITY MODEL FOR THE ELASTIC FORM FACTOR IN QUANTUM MECHANICS 

Factorization of hard form factors for large momentum transfers is the main relevant ingredient of the approach to the 
pion form factor advocated in Q . Hence, for any interaction of Coulomb-plus-confining type this model can be tested in 
quantum mechanics, which has already proven to be a rather efficient tool for studying various features of QCD (2l| - [25| . 

Within potential models, the elastic form factor of the ground state, F(Q), is given in terms of the wave function ^> by 

F(Q) = J d 3 r exp(iq • r) |*(r)| 2 = J d 3 fc*(fc) *(fc + q), Q = \q\, (4.1) 

where the bound-state wave function ^ is computed by solving the stationary Schrodinger equation for the Hamiltonian 

k 2 a 

H= — -- + V coni (r), r=\r\. (4.2) 
Ira r 

At large values of Q, the asymptotic behaviour of the elastic form factor F{Q) is given by the factorization theorem (26[: 



HO) — - F. K/) = 167r y^ 9 , R g ee |*(r = 0)| 2 . (4.3) 



The quantum-mechanical LD sum rules for decay constant R g and form factor F(Q) are rather similar to those in QCD: 



R g = I dkp^{k), (4.4) 

rfcrff(Q) /-feoff (Q) 

(4.5) 



F g {Q)R g = / / dfcidfeA2£(*i,fe2,Q). 

Jo Jo 



The spectral densities, p?^ and A^., are calculated from two- and three-point diagrams of nonrelativistic field theory; 
the derivation of these expressions up to 0(a s ) accuracy is presented in Appendix A. 



5 A way out would be to assume different effective thresholds for the O(l) and the 0(a s ) contributions to the pion form factor. This seems, 
however, a rather artificial construction. 
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The factorization limit (|4.3I) means for the momentum-dependent effective threshold in the three-point sum rule (|4.5[) 

fccff (Q) > k hD = (6tt 2 Rg) 1 / 3 . (4.6) 



For intermediate Q, the behaviour of F(Q) is controlled by the details of the confining interaction via the corresponding 
wave function ty. The LD model assumes that, also for intermediate Q, one may find a reasonable estimate for the form 
factor by setting fc e g (Q) — Hence, similar to QCD the only property of the bound state which determines the form 
factor in the LD framework is R g . To probe the accuracy of this approximation, we consider a set of confining potentials 



V con f(r) = a n (mr) T 



n = 2,l,l/2. 



(4.7) 



We adopt parameter values appropriate for hadron physics, i.e., m = 0.175 GeV for the reduced constituent light-quark 
mass and a = 0.3, and adapt the strengths <r„ in our confining interactions such that the Schrodinger equation yields for 
each potential the same *(r = 0) = 0.078 GeV 3 / 2 , which holds for a 2 = 0.71 GeV, o x = 0.96 GeV, and a 1/2 = 1.4 GeV. 

Figure [5] summarizes our findings in quantum mechanics. In the region Q < 2 GeV, the exact k e g(Q) is a complicated 
function of Q (Fig.[5]b). Moreover, its behaviour exhibits features similar to the exact QCD threshold, recalculated from 
pion form-factor data (see Fig. [3]). In the region Q > 2 GeV, the exact k e s (Q) is smoothly approaching its LD limit A;ld- 
For a steeply rising confining potential (n = 2), the LD model works with, e.g., 5% accuracy already for Q > 2 GeV. For 
a slowly rising confining potential (n = 1/2), such high accuracy is reached not before Q ps 6 GeV. However, we identify 
an important universal feature that does not depend on any details of the confining interaction: The accuracy of the LD 
approximation for the effective threshold, k e s(Q) ~ as well as the accuracy of the corresponding elastic form factor 
increase with Q in the region Q 2 > 4 GeV 2 . Figure [5] illustrates this observation for several confining potentials with 
different large-r behaviour. On the basis of these findings we are forced to conclude that, if at Q 2 ps 4-8 GeV 2 the LD 
model provides a good description of the data, the accuracy of this model won't become worse at larger values of Q 2 . 



Q 4 F(Q) [GeV 4 




Q [GeV] 



ketf(Q) [GeV] 
0.9 




Q [GeV] 



(b) 



Fig. 5: (a) Exact vs. LD ground-state form factors F(Q). (b) Exact vs. LD effective thresholds k c g(Q). Dotted (magenta) lines: 
local-duality limit; short-dashed (black) lines: harmonic-oscillator confinement, n — 2; full (red) lines: linear confinement, n = 1; 
long-dashed (blue) lines: square-root confinement, n = 1/2. 



The crucial question is to which extent these findings can be carried over to predictions for the pion elastic form factor 
by LD sum rules in QCD. In spite of subtleties related to evident differences between the spectral densities in relativistic 
and nonrelativistic theories [27j, we are convinced that the pivotal lesson from quantum mechanics also applies to QCD: 
Above the relatively low value Q 2 rs 4-6 GeV 2 , the accuracy of the approximation s e g (Q 2 ) ~ sld will increase with Q 2 . 
The results for the form factor found in clearly contradict this observation: the equivalent effective threshold for 

these models deviates from the LD threshold sld = 47r 2 f 2 in the broad range Q 2 « 4-10 GeV 2 and for QHj| even up to 
Q 2 pa 100 GeV 2 ; more importantly, the deviations between the equivalent thresholds and the LD threshold rise with Q 2 . 



5. CONCLUSIONS AND OUTLOOK 



We investigated the expected accuracy of the pion elastic form factor obtained by the LD version of QCD sum rules. 
The LD sum rule for the elastic form factor of some bound state may be constructed in any theory where hard exclusive 
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amplitudes satisfy the factorization theorem (that is, in essence, in any theory with interactions behaving Coulomb- like 
at small distances and confining at large distances). In this approach, the form factor is determined by two ingredients: 
the double spectral densities of the diagrams of perturbation theory and the Q 2 -dependent effective threshold s e g(Q 2 ). 
The effective threshold satisfies rigorous constraints at Q 2 = and for Q 2 — > oo but is unknown for intermediate values 
of Q 2 . The LD model assumes that approximating s e g (Q 2 ) by its value at Q 2 — > oo, 4ir 2 f 2 , yields reasonable estimates 
for the form factor at not too small Q 2 . We have tested the accuracy of the LD model in quantum mechanics, where the 
exact effective threshold may be calculated from the exact form factor, obtained by solving the Schrodinger equation. 
The new results reported in this work are the following: 

1. For Q 2 < 4 GeV 2 , the exact effective threshold exhibits a rapid variation with Q 2 , implying that in this region the 
accuracy of the LD model depends on the details of the confining interactions and cannot be predicted in advance. 

2. For Q 2 > 4-6 GeV 2 , independently of the details of the confining interactions, the maximal deviation of the exact 
effective threshold from the LD approximation is reached for Q 2 « 4-6 GeV 2 . As Q 2 increases further, the exact 
effective threshold approaches the LD limit quickly, thus improving the accuracy of the LD model for the elastic 
form factor. For generic confining interactions, the LD approach gives very precise results for Q 2 > 20-30 GeV 2 . 
In other words, it is not possible to predict at which value Q 2 of the momentum transfer the LD model provides a 
good approximation (with, say, 5% accuracy) to the exact form factor, as the precise Q 2 depends on subtle details 
of the confining interactions. However, one important conclusion may be drawn from our quantum-mechanical 
analysis: independently of any subtle details, the accuracy of the LD approximation increases for Q 2 > 4-6 GeV 2 . 

3. Existing data on the pion elastic form factor indicate that the LD limit sld = 4-7T 2 / 2 of the effective threshold may 
be reached already at the relatively low values Q 2 « 5-6 GeV 2 of the momentum transfer. 

On the basis of the results obtained in this study, we are forced to conclude that those large deviations from the LD limit 
in the range Q 2 = 20-50 GeV 2 reported in [J[| appear to us highly unlikely. 

Of course, our analysis does not provide a proof of but an argument for the accuracy of the LD approximation in QCD 
and it gives some hint towards the behaviour of the pion elastic form factor Ft T (Q 2 ) to be expected for large values of Q 2 . 
Accordingly, the accurate experimental determination of F^(Q 2 ) for Q 2 = 4-10 GeV 2 expected to be achieved by JLab 
will have important implications for the predicted behaviour of F 7r (Q 2 ) at large Q 2 up to asymptotically high momenta. 
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Appendix A: Perturbative expansion of Green functions in quantum mechanics 

We construct the perturbative expansions of both polarization operator and vertex function in quantum mechanics. 

1. Polarization operator 

The polarization operator H(E) is defined by [23| 

11(E) = (r' = 0\G(E)\r = 0), (A.l) 
where G(E) is the full Green function, i.e., G(E) — (H — E)" 1 , defined by the model Hamiltonian under consideration 

k 2 

H = H + V(r), H Q = —, r = \r\. (A.2) 

Am 

The expansion of the full Green function G(E) in powers of the interaction potential V has the well-known form 

G(E) = G (E) - G (E)VG (E) + G (E)VG (E)VG (E) + ■■■ , (A.3) 
with Gq(E) = (Hq — E) . It generates the corresponding expansion of H(E): 

11(E) =Il (E) + III (E) + --- . (A.4) 



Explicitly, one finds 



n o (£0 



ni(-B) 



{2nf £L-E' 

v ' 2m 



1 



d 3 fc d 3 fc' 



(2tt) 6 J £L-E¥r-E 



V((k-k') 



l\2\ 



2m 2m 

We consider interaction potentials V(r) which consist of a Coulombic and a confining part: 

V(r)=--+V amt (r). 



(A.5) 
(A.6) 

(A.7) 



Then the expansion (|A.4I) becomes a double expansion in powers of the Coulomb coupling a and the confining potential 
V con { (see Fig.©. 




+ 




E E 



+ 




n (E) Hj (e) 

Fig. 6: Expansion of the polarization operator in terms of Coulomb (wavy line) and confining (dashed line) interaction potentials. 

The contribution to Tl(E) arising from the Coulombic potential is referred to as the perturbative contribution, n pcrt . 
The contributions involving the confining potential V con f (including the mixed terms receiving contributions from both 
confining and Coulomb parts) are referred to as the power corrections, n powor . For instance, the first-order perturbative 
contribution reads 



u[ a) {E) 



1 



d 3 fc d 3 fc' 



4na 



(2n)*J »-Eg.-E{k-Wy 



am 
8^ 



d 3 k 



(A. 



The integral diverges but becomes convergent after applying the Borel transformation l/(a — E) — » exp(— aT). e The 
Borel-transformed polarization operator H(T) has the form 



n(T) = n pcrt (T) + n power (t) , n port (t) 



3/2 



2-kT) 



To? + 0(a 3 ) 

3 



(A.9) 



In the LD limit, that is, for T — >• 0, only Il pcrt (T) will be relevant. Nevertheless, as an illustration we provide also the 
result for the power corrections n power (T) for the case of a harmonic-oscillator confining potential V con f(r) = mui 2 r 2 /2: 



'<^(S 



3/2 



--w 2 T 2 | 1 



11 
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V 2-KniTa 



19 

480" 



w 4 T 4 



(A.10) 



Let us point out that n powor (T = 0) vanishes, similar to QCD. The radiative corrections in n pcrt (T') have a less singular 
behaviour compared to the free Green function, so the system behaves as quasi-free system. In QCD such a behaviour, 
frequently regarded as an indication of asymptotic freedom, occurs due to the running of the strong coupling a s and its 
vanishing at small distances. Interestingly, in the nonrelativistic potential model this feature is built-in automatically. 

Now, according to the standard procedures of the method of sum rules, the dual correlator is obtained by applying a 
low-energy cut at some threshold fc c ff in the spectral representation for the perturbative contribution to the correlator: 



1 

2^2 



dfc k exp 



2m 



-T 



irma (irmaY 



k 



3fc 2 



0(a 3 



IL 



c(T)- 



(A.11) 



Note that the Borel transform of the Green function (H — E) 1 yields the quantum- mechanical time-evolution operator in imaginary time 
U(T) = exp(-HT). 
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By construction, the dual correlator IlduaKT 1 , fcoff) is related to the ground-state contribution by 

n d uai(T,fc cff ) = n g (T) = i? g exp(-£ g T), R s = |V g (r = 0)| 2 . 



(A.12) 



As we have shown in our previous studies of potential models, the effective continuum threshold defined according to 
(|A.12j) is a function of the Borel time parameter T. For T — 0, one finds 



am 

4-7T 



(A.13) 



2. Vertex function 



We now calculate the vertex function T(E, E', Q), defined by 

T(E, E', Q) = (r> = 0\G(E)J(q)G{E')\r = 0), Q= \q\ 



(A.14) 



where J(q) is the operator which adds a momentum q to the interacting constituent. The expansions (IA.3I) of the full 
Green functions G(E) and G(E') in powers of the interaction entail a corresponding expansion of T(E, E', Q), cf. Fig. [7] 



T(E, E', Q) = T (E, E', Q) + T^E, E', Q) 



(A.15) 






EE E' E E' E' 

r (E,E',q) r^(E,E',q) 
Fig. 7: Nonrelativistic Feynman diagrams representing the lowest perturbative contributions to the vertex function F(_B, E' , Q). 



For the vertex functions Ti{E, E', Q), i = 0, 1, . . . , in (|A.15[) . their double spectral representations may be written as 

w z dZ ' l A l (z,z',Q). (A.16) 



T i (E,E',Q) = J 



The vertex functions Ti(E, E', Q = 0) and the polarization operators Iii{E) satisfy the Ward identities 

which are equivalent to the following relations between the corresponding spectral densities: 

lim Ai(z, z', Q) = S(z - z') Pi{z). 

Q->0 



(A.17) 



(A.18) 



We shall c onside r the double Borel transform E T and E' -> V: (a-E)- 1 -> exp(-oT), (a'-E')- 1 -> exp(-a'T'). 
Equation (|A.18I) leads to the following Ward identities for the Borel images: 



r i (T,T',Q = Q)=U i (T + T'). 



(A.19) 



a. One-loop contribution To to the vertex function 
The zero-order one-loop term has the form 

T (E,E',Q) 



1 



d 3 fc 



k- 



E 



( (k+q) 2 
\ 2m 



- E' 



(A.20) 
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This may be written as the double spectral representation 

T (E,E',Q) -- 

where 



_* E -2- E 

2m ^ 2m ^ 



7 A (z,z',Q), 



(A.21) 



A (z, z', Q) = ^ j d 3 k6(z - k 2 ) 5(z' -(k-q) 



1 7T 



(27r)32Q ^'--Q 2 ) 2 - 4^<0). (A.22) 

Hereafter, we use the notations k = y/z and k' = \fz' '. In terms of the variables k and k', the 9 function takes the form 

6((z' -z- Q 2 ) 2 - 4zQ 2 < 0) = 9(\k -Q\<k' <k + Q). (A.23) 



b. Two-loop contribution Ti to the vertex function 



We consider here only corrections related to the Coulomb potential, since power corrections induced by the confining 
interaction vanish in the LD limit. The two-loop 0(a) correction receives two contributions and has the form 



1 



d 3 fc d 3 fc' 



47ra 



Vl{E ' E '- Q) WJ ^-E^-E'(k-(k'-q)) 2 



1 



1 



(fc'-q) 2 



E> 



(fc+q) 2 



E 



(A.24) 



Having in mind the subsequent application of a double Borel transformation in E and E' , it is convenient to represent 

with 

d 3 fc' f d 3 k 



Ti as a sum of two terms, T\ = 1^ + with 



r[ a \E,E',Q) = ^ 



am 
8^ 
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(A.25) 



The double Borel transformation in E — > T and _E' — > T' is now easily performed. 
Let us start with . One integration in may be performed, leading to 



T[ a \E,E',Q) 



16tt 2 



d 3 fc' 



d 3 k 



(k-q) 2 



, -E)(f--E')\k\ 

2m / V 2m / 1 1 



(A.26) 



The first term corresponds to the contribution of the "left" two-loop diagram in Fig. [71 i.e., with the potential before the 
interaction with the current J(q), while the second term is represented by the "right" two-loop diagram in Fig.[7J The 
corresponding double spectral densities have a form very similar to Ao: 

A# (k, k\ Q) = 0(\k - Q\ < k' < k + Q) 9(0 < k) 9(0 < k'), A^(fc, k', Q) = A&> (k\ k, Q). (A.27) 

Explicit calculations yield the following double spectral densities of the two contributions to 1^ related to the "left" 
and "right" two-loop diagrams in Fig. [7j 



Afl(k,k',Q) = 



am 1 
32tt 6 Qk 



log 2 



k' -Q + k 



k' -Q-k 



-log 2 



k' + Q + k 



k' + Q-k 



A^(k,k',Q) = A\°>(k',k,Q). (A.28) 



At Q = 0, r[ a) (T,T',Q = 0) satisfies the Ward identity, r[ a) (T,T',Q = 0) = n[ a) (T + T'), whereas r[ b) (T,T' ,Q = 0) 
vanishes: ) (T, T', Q = 0) = 0. For large Q and T, T' ^ 0, Tf (T, T' , Q) assumes a factorizable form (see Eq. (|A~25)) ): 



rf(T,T',Q)^^=n (T)n (TO. 



(A.29) 



At the same time, both Tq(T, T', Q) and (T, T', Q) are exponentially suppressed for large Q and T, T' ^ 0. Hence, 
determines the large-Q behaviour of the vertex function. 
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Dual correlator 



The dual correlator Td-aal(T, T", Q) is constructed in a standard way, by application of a low-energy cut to the double 
spectral representation of the perturbative contribution to (|A.16|) : 



k eH (Q,T) fceff(Q,T') 

u2 \ r / f„'2 



r d „ai(T,T',Q)= y* dfc2fcexp(-^-r) jf dfc' 2fc' exp(-^- T') A(z, z', Q) + r powor (T, T', Q). (A.30) 







By construction, the dual correlator corresponds to the ground-state contribution exp(— E g T) exp(— E g T') R g F g {Q). 

In the LD limit T = and T" = 0, r powcr (T, T", Q) vanishes and the ground-state form factor F g (Q) is related to the 
low-energy part of the perturbative contribution considered above: 

fccff(Q) fceff(Q) 

J dk2k J dk'2k'A(k,k',Q) = F g (Q)R g , (A.31) 
o o 

with A(fc, k', Q) = A (fc, k', Q) + A$ (k, k', Q) + A^ (k, k', Q) + A^ (fc, k', Q) + Ag> (fc, fc', Q) + 0(a 2 ). 

In order to provide the correct normalization F g (Q = 0) = 1 of the elastic form factor F g (Q), the effective thresholds 
should be related to each other according to k e s {Q = 0) = fc ff ; then the form factor is correctly normalized due to the 
Ward identity (|A.18|) satisfied by the spectral densities. 
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